Time-Local Quantum-State-Diffusion Equation for Multilevel Quantum Dynamics 
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An open quantum system with multiple levels coupled to a bosonic environment at zero tempera- 
ture is investigated systematically using the non-Markovian quantum state diffusion (QSD) method 
[W. T. Strunz, L. Diosi and N. Gisin, Phys. Rev. Lett. 82, 1801 (1999)]. We have established exact 
time-local QSD equations for a set of interesting multilevel open systems including the high-spin 
systems, multiple transition atom models, and multilevel atom models driven by time-dependent 
external fields. These exact QSD equations have paved a way to evaluate the dynamics of the open 
multilevel atomic systems in the general non-Markovian regimes without any approximation. 
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I. INTRODUCTION 

Many quantum processes, including emission, absorp- 
tion and quantum interference [l|, 0|, typically involve 
multilevel atomic systems coupled to a quantized bosonic 
field. For example, quantum multilevel models are used 
in studies of the high-spin model [3j, super-radiance, 
quantum phase transition in Dicke model 0. Q , electro- 
magnetically induced transparency (EIT) [6(, molecular 
aggregate |7[ and many other phenomena in quantum 
optics [1, Hi and quantum chemistry [Hldlj]. More re- 
cently, multilevel atomic systems also play a crucial role 
in describing coherent state transfer [12| and entangle- 
ment dynamics and control [l3l - [T8j . All of those models 
concern the interactions between a multilevel atomic sys- 
tem and its surrounding environment and external driv- 
ing fields leading to collective and competitive behav- 
iors including dissipation, fluctuation, decoherence and 
revival. Approximate master equations and other alter- 
natives such as the Langevin equations or quantum tra- 
jectories are used to describe quantum dynamics of the 
multilevel atomic system when the weak-coupling and 
Markov approximations are valid [l9j . 

In this paper, the multilevel atomic systems will be 
treated in the framework of the non-Markovian quantum 
state diffusion equation [20|-[25j without any approxima- 
tions. Unlike the methods of Markov quantum state dif- 
fusion or quantum jump simulations [26l430l |. where quan- 
tum trajectories are unravelings of the density operator 
of the system of interest, the non-Markovian QSD equa- 
tion is derived from the first principle and is determined 
uniquely by the system Hamiltonian, the coupling op- 
erator (called the Lindblad operator) and the spectral 
density of the environment [2l|, [3l[ . It has been shown 
that the non-Markovian QSD equation can be useful in 



the following ways: (i) it is an exact description of quan- 
tum dynamics of the system under the influence of the 
environmental noise with finite memory time (32rl36| . in 
particular, it serves as a powerful tool in numerical sim- 
ulations of quantum open systems; (ii) it could also be 
used to derive the exact or approximate master equations 
for the reduced density matrix [37l - l39| . 

Exact QSD equations have been found in several inter- 
esting models including a two-level atom in the dissipa- 
tive environment [22| , a harmonic oscillator in Brownian 
motion [38j . a three- level atom in the dissipation model 
[25| . and a two-qubit model [4(|. However, it remains 
unknown how exact QSD equations can be derived for 
general multilevel systems coupled to a dissipative en- 
vironment. The purpose of this paper is to establish a 
set of exact QSD equations for a large class of multilevel 
atomic models coupled to a dissipative environment with 
or without external driving fields. 

This paper will be organized as follows. In Sec. HH 
we will begin by reviewing the basic concepts of the 
non-Markovian QSD approach that will be useful for the 
discussions presented in the subsequent sections includ- 
ing the formal QSD equation and the O-operators. In 
Sec. IIII1 we will then proceed to discuss the construction 
of the time-local QSD equation by explicitly determining 
O-operators in various multilevel models. Afterwards, 
we present the numerical results from QSD simulations 
to illustrate the non-Markovian dynamics of multilevel 
atomic systems under environmental noises in Sec. IIVI 
Finally we will conclude the paper in Sec. |V] 



II. NON-MARKOVIAN QSD EQUATION 

In the framework of system-plus-environment, the to- 
tal Hamiltonian is given by (setting h — I): H to t = 



H sys + Ek(ffk La k + Pk^Ok) + Ek Uk a t a k, 



where H sys is 
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the system Hamiltonian and L is said to be the Lindblad 
operator coupling the system to the environment. The 
environment is described by a set of harmonic oscillators 
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a^ak satisfying [akjfl£/] = <5k.k'- The non-Markovian 
QSD approach is designed to recover the reduced density 
matrix for the open system by the average of the pure 
states driven by a colored Gaussian process z£: 

Pt = M[\Uz*))(Uz*)\l (1) 

where M stands for the statistical average over the noise 
z t *. The dynamics of the stochastic unravelings ipt{z*) 
(quantum trajectories) is governed by the non-Markovian 
QSD equation [UIH|]: 

rt 



-iH, 



sys 



Lz t - ^ f / dsa(t,s)— ip t (z*), 

JO 0Z s . 

(2) 

where a(t,s) = ^ k |gk| 2 e 1 * s ' is the environmental 
correlation function. By construction, the Gaussian noise 
z t * satisfies M[z$] — M[z^z*] = and M[z t z*} = a(t,s). 
Note that the functional derivative contained in Eq. @ 
is a time non-local term depending on the entire evolu- 
tion history from to t. This time non-local term is a 
major obstacle in using Eq. J2} as a numerical tool or 
as an analytical approach to deriving the corresponding 
non-Markovian master equation. One way to transform 
the formal non-Markovian QSD equation @ to a time- 
local equation is to introduce an O-operator satisfying 
H *(f } = 0{t, s, z*)tp t {z*), then the QSD equation could 
be recast into a convolutionless form: 



-iH sys + Lz* t ~ tfQ(t, z*)} il> t {z*), (3) 



where 0(t, z*) = f Q dsa(t,s)0(t,s,z*) and the O- 
operator may be determined from the following equation: 

dO(t,s,z*) 



dt 



-iH sys + Lz* t -L*0(t,z*), 



0(t,s,z*)} - L 



^50(t,z* 



Sz* 



(4) 



Hence, the key issue in establishing the exact time- 
local QSD equations is to explicitly construct the O- 
operator defined in (j4j. In what follows, we shall show 
how to determine the O-operator explicitly for a large 
class of multilevel atomic systems. 

It is worthwhile to note that, from the density ma- 
trix reconstruction defined in Eq. (JTJ and Novikov theo- 
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371 ] . we may obtain an exact master equation from 
Eq. ©: 

d tPt = [-iH sys , pt] + [L, M[P t d*}} + [M[OP t },L% (5) 

where P t = \tpt(z*))(ipt(z*)\. The master equation takes 
a more concise form when the O-operator is noise-free: 
0(t,s,z*) = 0(t,s), 

d tPt = -i[H sys , Pt ] + [L lPt O f ] + [Opt,L\ (6) 

In particular, under the Born-Markov approxima- 
tion with the delta function a(t, s) = T6(t — s), then 
0(t, s, z*) — TL/2. Eq. §6§ recovers the Lindblad master 
equation: 

d t pt = -i[H Bys ,p t ] + ^ ([L,pttf] + [Lpt,tf]) . (7) 



III. TIME-LOCAL QSD EQUATIONS 
A. High-Spin Model 

Now we consider a collective angular momentum 
(Dicke) model or the high-spin model described by 



u>J z . L = J_, 



(8) 



where J z and J_ are the angular momentum operators 
with spin-L It can be proved (see Appendix [5} that 
the O-operator contains 1(21 + 1) terms with the basis 



operators (ffl = \ j) (j + k + 1| , j = 1 , 



0. 



,21 — 1. And the terms with O 



,21 

(k) 



-k, for each 
contain fcth 



order noises for k = 0, • • • , 21 — 1. 

Clearly, one may construct many different sets of basis 
operators, but they can always be realized by linear com- 



binations of the operators O 



(fc) 
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operators may also be given by J 



selected here. The basis 

: 1 ' : . The 
simplest model is the spin-1/2 system (I = 1/2), where 
the O-operator contains only one term with the basis op- 
erator = er_ . Another model is spin-3/2. We can 
show the O-operator 0(t, s, z*) totally contains 6 terms, 
amongst, only one term contains the second-order noise: 



O 



3 

rt 



fj(t,s)0 



(0) 
3 



2 

E 



(2) 
Pi 



(t, s, si,s 2 )z* Sl z* S2 dsids 2 0\ 



(2) 



(9) 



In Table HI we list the numbers of the basis operators 
for the high-spin models from spin-1/2 to spin-7/2. 
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TABLE I. The number of basis operators with fc-fold inte- 
gration over noises in the O-operator for spin-Z system in a 
dissipation model. N is total number of the terms in the 
O-operators. 

The Hamiltonian and Lindblad operator considered in 
Eq. © imply that the energy levels of the multilevel atom 
are equidistant and all the exited levels dissipate to the 
lower neighboring levels. However, these two constraints 
can be relaxed. We start with the model Eq. ([5]) with 
the modified Hamiltonian and coupling operator: 



H sys =yiojj\j){j\, 
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L = J2^\j)(j + l\ 



(10) 



More specifically, we use a three-level system as an exam- 
ple for simplicity. It can be verified that the O-operator 
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can be constructed as follows (For details, see Appendix 
H]or Eq. (6) of H3): 

= f 1 (t,s)J_+f 2 (t,s)J z J^ + f p(t,s,s 1 )z* 1 ds 1 Jl. 

Jo 

(11) 

Consequently, 0(t,z*) = Fi(t)J_ + F 2 {t)J z J- + 
JlP{t,s 1 )z* 1 ds 1 J'i, where Fj(t) = f*a(t i a)f j (t,s)ds, 
j — 1,2 and P(t,si) = f* a(t, s)p(t, s, Si)ds. The ini- 
tial conditions for these coefficients are fi(s, s) = k 2 /V2, 
f 2 (s,s) = (n 2 — Ki)/V2, p(s,s,sx) — and they satisfy: 

d t fi = i(oj 3 - u 2 )fi - \/2[(«i - «2)-Fi - wi-Fbl/x 

- V2iKiP(t, s), 

dth = i(^3 + wi - w 2 )/i - V2iHnP(t, s) 
+ i{uj 2 - wi)/ 2 + V2/si(.Fi - F 2 )/ a 

- v^[(2/si - k 2 )F x - 2 Kl F 2 ]h, 

d t p = i(ui 3 - u)i)p + V2n 2 Fip + V2kiP(/i - f 2 ), 

p(t, s, t) = i^^f^t, s) - *^/ 2 (t, s). (12) 

These equations will reduce to Eqs. (7), (8) and (9) in (2f| 
by setting W3 = oj, lu 2 = 0, ui% = —uj and ki = k 2 = \2. 
As shown in the next subsection, the time-local QSD 
approach can be extended to the driven atomic models 
with the multiple transitions between the energy levels. 

B. Multiple Transition and Driven Atomic Models 




FIG. 1. (Color online) Schematic diagram of a general EIT 
model: A four-level atomic system with two driving fields 
(green double arrow lines) as /S.2(t)T23 + h.c. and A4(t)T34 + 
h.c The admitted damping channels (blue single arrow lines) 
are T 12 , T 13 and T M . 

In this subsection, we consider a driven multi-level 
atom with multiple dissipative channels. The system 
Hamiltonian is given by, 

N 

H sys = J2"j\j)(j\+ E [&k lk2 (t)\h)(k 2 \ + h.c.], (13) 

3 = 1 k t ^k 2 



where Afc 1 fc 2 (t)'s are the time-dependent functions and 
1 ^ ki,k 2 ^ N. Also, we consider a general Lindblad 
operator, which is given by 

L= 2J Knxn 2 \n\){n 2 \. (14) 

We remark that one may always determine the exact 
time-local QSD equation for the above generic Hamil- 
tonian and the Lindblad operator. However, it can be 
verified that a particularly simple noise-free O-operator 
exists if the following two conditions are satisfied: (i) 
any operator |fci)(fc 2 | appearing in the driving term is 
not present in the Lindblad operator in Eq. (fT4|) : (ii) 
the cycle transition terms, for example, |ni)(n 2 |, ^2) (^3! 
and |m)(n3|, cannot be simultaneously contained in the 
Lindblad operator. As an illustration, we consider an in- 
teresting model consisting of a four-level atom shown in 
Fig. ([1]), the system is described by 

4 

H sys = 5>ili>01 + [A 2 (*)|2)<3| + A 4 (£)|3>(4| + h.c], 

L = k 2 |1)(2| + k 3 |1)(3| + k 4 |1)(4|, (15) 

where A 2 (t) and A 4 (i) are independent driving exter- 
nal fields. Thus the O-operator can be constructed as 

0(t,s,z*) = E}=a/j(M)|l>0'l with fj( s > s ) = k j- And 
the differential equations for the coefficients are 

4 

dth(t, s) = -iwia/2 + i&lh + K *jfi> 

J'=2 

4 

dth(t, s) = -iu, 13 f 3 + iA 2 / 2 + iA*J 4 + F 3 J2 rffj 

3=2 

4 

dth{t, s) = -iwu/i + iA 4 / 3 + F 4 J2 K jfi> ( 16 ) 

where u)\j = U\ —cjj and Fj = Fj (t) = J* dsa(t, s)fj (t, s), 
J = 2,3,4. 

When the driving terms in H sys are omitted, the model 
in Eq. (| 13[) reduces to an JV-level atom with multiple 
transition channels (4lj | where the transition takes place 
between the highest level and all the lower energy levels, 
and transitions between any other levels are forbidden. 
Explicitly, the Hamiltonian and Lindblad operator in this 
case are given by 

N JV-l 

H sya = J2"i\3)(j\, L=J2^\j)(N\, (17) 

respectively. For this model, we can show that 
the O-operator can be explicitly constructed as 
0(t,s,z*) = E^i 1 fj(t, s)\j){N\. Hence 0(t,z*) = 
^2^=1 Fj(t)\j}(N\ with the initial conditions, fj(s,s) — 
Kj , j = 1, • • • , JV — 1. By the consistency condition Q , 



4 



one gets 

N-l 

dtf.it, s ) = i(u N - Ui )fi + fi £ ****(*)• (is) 

fc=i 

In a more general case, we may consider an [N + M)- 
level atom with an upper energy band consisting of N 
levels and a lower band consisting of the other M levels. 
We assume that transitions between the upper band and 
the lower band are allowed, but those between the energy 
levels inside the upper band or the lower band are strictly 
forbidden. Such a model may be described by 

N+M M N 

H sys = "MUl L = J2 E Kjk\j)(k\. (19) 

j=l j=l k=M+l 

The O-operator is explicitly constructed as a noise-free 
formation by 0(t, s) = £V =1 Y,k=M+i fjk(t, s)\j)(k\ and 

6(t) = J2f=lJ2k=M+l F jk(t)\j)(k\ With f jk (s,s) = K jk . 
And we have 

M N 

9tfjk{t, s) = i(uj k - Uj)f jk + E E fjk'n}k' F j'k- 

j'=l k '=M+l 

(20) 

Therefore, for the cases with the noise-free O-operators, 
the exact master equations can be derived directly from 
Eq. ([5]). It is noted that such master equations may not 
be of a standard Lindblad form, but their positivity is 
automatically guaranteed by the derivation. 

IV. NUMERICAL RESULTS AND 
DISCUSSIONS 




FIG. 2. (Color online) Coherence and von Neumann entropy 
[4^ 1 of a dissipative four-level model: H sys = U)J Z , and L — 
J_, and I = 3/2. Initially, Vo = (1/2)(|1> + |2) + |3) + |4)). 
The time evolution of coherence | {pi4 ) | and entropy S v n = 
—p log 2 p is averaged over 1000 trajectories for different values 
of 7 (r = id) with the first-order noise term. 

It is known that the stochastic wave- function ip t (z*) in 
Eq. ([3]) does not conserve the norm of the wave-function. 
To efficiently simulate quantum open systems, one usu- 
ally uses the nonlinear QSD equation for the normalized 



state Vt = n^jy: 

jfo = [ - iH sya + (L— (L) t )r t + (l% (o(t, r ) 
- (0(t, s*)) t ) - (ito( tj r) - (it<5(t, r)) t ) 

(21) 

where (L) t = (tjj t \L^\'ijjt} and z* t = z t * + f* a*(t, s)(L j <) s ds 
is the shift complex Gaussian process. 

The non-Markovian QSD approach is valid for an ar- 
bitrary correlation function. For simplicity and consid- 
ering Markov limit, the non-Markovian effect is modeled 
by the Lorentz spectral density: S(ui) = 2 r ^ 2 . Then 
the correlation function obtained from the Fourier trans- 
formation is given by a(t, s) — / °° dijjS{ijj)e^ lu) ^ s ^ = 
_JZ e -7|t-s| ; where I/7 is an important non-Markovian 
parameter representing the memory time of the environ- 
ment. When 7 — > 00, the correlation function approaches 
the Markov limit with a(t, s) —> TS(t — s). 

As our first example, we consider the numerical sim- 
ulation of the dissipative dynamics of a four-level sys- 
tem with the O-operator given by Eq. (0). In the case 
of dissipative bath at the zero temperature, the quan- 
tum coherence of the four-level system will decay and 
the populations of the excited levels will be transferred 
to the ground state with time. When the bath is in a non- 
Markovian regime with 7 = 0.3, it is shown that the deco- 
herence dynamics deviates from the exponential decay of 
the Markov case (See the bump of pu in Fig.[2fa)). Con- 
sequently, the decoherence is delayed as the entanglement 
between system and bath builds up. We also plot the von 
Neumann entropy [5 v n = — Tr^lnp)] (setting the Boltz- 
mann constant ks = 1) in Fig. HJb). When the parame- 
ter 7 increases to 1.0 (moderate non-Markovian regime) 
and 3.0 (near-Markov regime), the coherence quickly de- 
cays to zero at the time point (uit ~ 5). However, we 
see that their entropy curves are different. For example, 
when 7 = 3.0, SVn approaches zero at tot ~ 5, which 
means that the decoherence time between the ground 
state and the highest energy level and purification time 
(That is, iSvn — 0) coincide. Interestingly, we can see the 
5 v n is significantly modified by the non-Markovian effect 
with 7^1. 

In Fig. [31 we show the dynamics of a driven four-level 
atomic system coupled to a dissipative environment with 
different 7's. When t — 0, the atom is totally populated 
at the highest level |3). With the modulation of the driv- 
ing fields A2(f) and A^t), Level 2 and Level 4 will be 
coupled to Level 3, which helps to establish the coherence 
terms P23 and P34. Those coherence terms will decay due 
to the dissipative channels to Level 1. From Fig. EJa) 
and Fig. [3fb) , when 7 = 0.3 the non-Markovian environ- 
ment with long memory time clearly increases the mag- 
nitude of the coherence terms. It should be noted that 
the differences between 7 = 1.0 and 7 = 3.0 are not sig- 
nificant since the system is relatively close to the Markov 
regimes. We can also see from Fig. [3Jc) and Fig. [2Jd), 
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FIG. 3. (Color online) Coherence and population dynam- 
ics of a dissipative four-level model driven by two external 
fields. The model is described by Fig. [T] and Eq. US]). The 
parameters are chosen as follows: u\ = O.fr, u)2 = 0.3r, 
uj'i = 0.6T, uj 4 = 0.2r, A 2 (t) = 0.lFe 2it T 2 3 + h.c, A 4 (t) = 
0.ire 2it T 34 + h.c, k,2 = 0.4, k 3 = 0.8 and Ka = 0.3. Initially, 
we choose tpo = |3). The time evolution of the reduced den- 
sity matrix is obtained by averaging over 1000 trajectories for 
each 7 obtained from the QSD equation with the first-order 
noise term: (i) 7 = 0.3, the blue solid lines; (ii) 7 = 1.0, the 
red dashed lines; (iii) 7 = 3.0, the black dot-dashed lines. 



that when 7 = 0.3, the non-Markovian effect causes pop- 
ulation ffuctuation of Level 3 and Level 1. Typically, we 
see that the population transfer rates to the ground state 
are increased by the shorter environment memory time. 



V. CONCLUSION 

In summary, we have studied the non-Markovian dy- 
namics of multilevel atomic systems using the non- 
Markovian quantum trajectory method. The time-local 
QSD equations are obtained by the explicit O-operator 
construction for several physically interesting models. As 
shown in this paper, the time-local O-operators allow nu- 
merical simulations to be implemented efficiently for mul- 
tilevel open systems. For the multilevel atoms without 
driving fields, we show explicitly how to construct the 
exact O-operators containing a finite-order noise terms. 
For the atomic models with multi-transition channels in- 
cluding the models with driving external fields, we verify 
that the O-operator can take a noise-free form, so the ex- 
act QSD equations can be easily established. In the case 
that the O-operator contains no noise, the exact convo- 
lutionless master equation can be derived. The results of 
this paper will open am avenue to exploring novel and fas- 
cinating phenomena in multilevel atom-field interaction 
problems where the non-Markovian features are impor- 
tant. 
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Appendix A: Proof for the O-operator in Dissipative 
High-Spin Models 

This appendix provides the details of deriving the 
O-operator for a multilevel atom dissipative model in 
Eq. ([8]). The techniques with some necessary yet not 
complicated modifications can be used to derive the O- 
operators for a general multilevel models including the 
models with the external driving fields. For the high- 
spin model, the system Hamiltonian and the Lindblad 
operator may be generally rewritten as 

21+1 

H sys = ]T C m H$ , = \m) (m\, 

m— 1 
2/+1 

L=J2 G nL n , L n = \n-l)(n\. (Al) 

71=2 

When C m 's and G„'s are chosen as the Clebsch-Gordan 
coefficients, these operators are referred to the angular 
momentum operators of spin-L In practice, Eq. (|Aip 
can be generalized to a general dissipative model for a 
multilevel atom with an arbitrary energy distribution. 
The model can also relaxed to the case containing time- 
dependent damping coefficients. 

In what follows, we will show that there are (21 — k) 
terms containing fcth-order noise (k = 0, • • • , 21 — 1) in 
the O-operator. Explicitly, 0(t, s, z*) may be expanded 
as 

21 21-1 r t 

nt ft 

H \- ■ ■ Pi (M>Sl: ' ' ' ,S2l-l)z* t ■ ■ ■ 

Jo Jo 

z* S2l _ 1 ds 1 ds 2 ---d S2l - 1 0f- 1 \ (A2) 

where the coefficients of (1 ^ j ^ 21 — k) are sym- 
metric functions of Si, ...Sfe. Subsequently, 

21 21-1 r t 

+ ■■■+[ ■■■ f p[ 2l ~ l) (t, ,»«_!)*:,••• 

Jo Jo 

z* S2i _ i d Sl ds 2 ---ds 2l _ 1 0? l - 1) . (A3) 
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>(*0 



In accord with definition, Fj(t) = J dsa(t, s)fj(t, s) and fj's and p^ s. When n in L n runs from 2 to 2/ + 1, 

Oj s will present themselves successively with j in 

(k) 

O) running from 1 to 21 — k. It is also a neces- 



• ,2/- 



jSfc) = f*dsa(t, s)pf\t, s,si,--- , s fe ), fe 



Clearly, the operators 0} form a set of basis opera- sary requirement for the initial condition of O-operator 



tors: 



O 



O) 



\j)(j + k + l\. (A4) 



Hence the total number of the basis operators O) is 
(21 + 1)1. In fact, any linear combinations of operators 
in Eq. (|A4I) may be used to construct the O-operator as 
long as they are linear- independent and satisfy Eq. (j4|). 
Here we use Eq. (|A4[) for the simplicity in the following 



proof. More often the notation O 



(fc) 



Jrv: 



1 jk+l 



has 



a more transparent meaning for the angular momentum 
model. 

By substituting Eq. (|A1J1 into Eq. g]), the left hand 
side of Eq. dU becomes 



2Z (-i 2/ 1 i- /, f <-\ 



pf'(t,s,s 1 )z* 1 dsi 



+ z* t pf\t,s,t) 

(2i-l)/. 

+ (22 - IK 



(i) 



(2i-l) 
Pi 



dt 

dsi--- ds 2 i-i 



(t,s,si 



,S2l-2,t) 



z* S2l _ 2 ds! 



ds 2 l-7 



(21-1) 



(A5) 



The right hand side of the equation consists of four terms. 
It can be easily seen that each of them can be expressed as 
a linear combination of the operators given in Eq. (|A4|) . 
The construction of the O-operator can be summarized 
in the following four crucial observations: 

(i) [— iH a y S , 0(t, s, z*)] consisting of [H^s, Oy ] can be 
decomposed as 

+ k + l|(<Sj.rn — Sj+k+l.rn) = 0^(Sj. m ~ Sj+k+l.rn)- 

(A6) 

(ii) [Lzj: ,0(t, s, z*)] consisting of the noise term 
z%[L n ,Oj] turns out to be 

k 



^(li-i)O' 



l\8j,n 



zKO^X'5^ - 0\ K ^S j+k , n _ 2 ). (A7) 



These terms correspond to those with z\ in Eq. (|A5I) . 
which will appear in the boundary conditions between 



0(s,s,z*)=L = J2 2 j+ 2 1 G n L n . 

(hi) [-tf6,0] is consisted by [of\\n)(n- l\of'\ 
which turns out to be 



U " .J))+k>* '~ U n'+l d j'+k'+l,j- 



(A8) 

(iv) For — iJ — g '* ; we only need to consider the 
terms with k ^ 1. The differential functional leaves the 
operator unchanged. Then one typical component is 

- Lt0f =-\j + + k + l\S n , j+1 = -Of^. (A9) 

It contributes to the presence of Oj s. 

To summarize the last four steps, we have proved the 
existence and constitutions of the O-operator for the an- 
gular momentum model of the multi-level atom. By in- 
troducing Eq. (|A4[) , the consistency condition equation 
(@| is shown to be consistent and complete. 

Now we consider a general four-level atom dissipative 
model with I = 3/2 in Eq. (jAl[) . The group of closed 
differential equations for the coefficients can be derived 
through the above analysis. They are given by, 



d t fi = 


-*(Ci 


- C 2 )h F 




9th = 


-*(C a 


- C 3 )f 2 4 


GH2F2 - G^/ 2 Fi - G* 2 p[ 1] 


dth = 


-i(C 3 


- C 4 )/ 3 4 


G4/3F3 — G* 3 f 3 F 2 — G* 3 P 2 ^ 


hp? = 


-i(Ci 


- cM 1] 


+ G* 2 hpW +G* 3 F 2 p£\ 


Vft w = 


-i(C 2 


C4)P 2 1] 


+ G* 3 f 2 pW + GtF 3 pV 



- G* 2 F lP2 1] - G* 2 f 3 PP - 2G* 2 p[ 2 \ 
d tP {2) = -t(Ci - C 4 )p1 2) + GlhP {2) + G* 3 pPpP 



+ g* 4 f 3P y 

Here the boundary conditions are, 



(A10) 



2pf\t,S,Si,t) = G 2 py(t,S,Sl) - GiP^\t,8,8l), 



p^ ) (t,s,t) = G 2 f 2 (t,s)-G 3 f 1 (t,s), 
p ( 2 1 \t,s,t) = G 3 f 3 (t, s )-G 4 f 2 (t,s). 



(All) 



The spin-3/2 dissipative model [See the O-operator in 
Eq. ©] can be solved by setting C m = — I — 1 + m, 
m = 1, 2, 3, 4 and G 2 = G 4 = V3, G 3 = 2. 
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